Body-rock or lift-off in flow 

Frank T. Smith*and Phillip L. Wilson^ 
January 14, 2013 



Abstract 

Conditions are investigated under which a body lying at rest or rocking 
on a solid horizontal surface can be removed from the surface by hydrody- 
namic forces or instead continues rocking. The investigation is motivated 
by recent observations on Martian dust movement as well as other small- 
and large-scale applications. The nonlinear theory of fluid-body inter- 
action here has unsteady motion of an inviscid fluid interacting with a 
moving thin body. Various shapes of body are addressed together with a 
range of initial conditions. The relevant parameter space is found to be 
subtle as evolution and shape play substantial roles coupled with scaled 
mass and gravity effects. Lift-off of the body from the surface generally 
cannot occur without fluid flow but it can occur either immediately or 
within a finite time once the fluid flow starts up: parameters for this are 
found and comparisons are made with Martian observations. 



1 Introduction 



The effect or the use of a flow of fluid to remove a body originally stationary or 
rocking on a fixed solid surface has many applications and is of interest over a 
range of length and time scales. The force required to do this washing, clearance 
or erosion of the original surface is also of much concern. 

The applications vary from removal of debris, erosion of soil by wind, water 
or raindrops, sand and pebble movement on beaches, dust loss, dust blowing, 
leaf-blowers and related geological and industrial phenomena, through cleaning 
and washing processes such as with fluid knives and cameras, conveyor de- 
sign, biomedical problems, weather damage, to aircraft take-off, estimation of 
runway length, sports applications such as ski-jumping, car and cycle racing, 
and the safety of wind-blown buildings. See for example Godone and Stanchi 



2011 , Bascom 1980 , Virmavirta ct al. 200f , Witt et al. 



f999 



In geolog- 
ical settings this process is known as saltation. The typical scales range from 
the comparatively small sizes of various biomedical settings, dust problems and 
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household/industrial cleaning to the larger sizes of tornados and tsunamis lifting 
relatively large obstacles Hunt 2005 Mikami et al. 2012 , and related disas- 



ters, along with military applications, and even on to planetary size. Moreover, 



Cox et al. 2012 , Hall et al. 2010 and references therein provide observational 



evidence that shoreline boulders whose weight is of the order of 10 tonnes can 
be moved by wave motion even in years without exceptional storm activity, 
and that boulders up to 78 tonnes have been moved to 11m above high wa- 
ter in significant storms, with the role of tsunamis ruled out. A comparison 
between boulder transport by storm waves and tsunamis is given in Barbano| 
et al. 2010 . In these geological applications, some of the more impressive feats 



of wave transport occur with long, slender boulders conforming to the assump- 
tions of the present study. Further, shallow granular avalanches — such as those 
occuring in pyroclastic flows and snow avalanches — are studied in Johnson and 



Gray 2011 , Gray and Ancey 2011 where attention is drawn to the effects of 



stationary granular material within evolving flows, in addition to pile collapse 
behaviour and segregation rates between constituents, where large and small 
particles percolate to the base or surface of an avalanche. In the human body, 
the lift-off due to hydrodynamic forces in the blood flow of thromboses formed 
near the site of blood-vessel injury, such as that caused by atherosclerosis, is 
of significant concern as the thrombus or its resulting transported embolus can 
lead to stroke, myocardial infarction, deep vein thrombosis, and other damaging 
conditions Ku 1997 . Hemodynamic lift-off and transport of tumour cells from 



damaged vascular walls may be one of several rate-limiting factors in cancer 

A dramatic example from human engi- 



Koumoutsakos et al. 2013 



metastasis 

neering and exploration was the loss of the space shuttle Columbia as a result of 
its collision with a relatively small piece of insulating foam separated from a fuel 
tank apparently by the flow past the shuttle Gomez III et al. 2003 . We should 



also mention human experience and pleasure, for example in blowing spilt salt 
off a table, and watching "dust devils" dance across a beach or even Mars. 

The issues involved include the need for understanding of the influences of 
length and time scales, the effects of the incident flow direction and velocity 
as well as the effects of the detailed obstacle shape, and issues of stall, lift, 
impact and drag. The roles of the Froude or Richardson, Reynolds and effective 
Stokes or Womersley numbers in particular need consideration. A related issue 
is the question of whether or under what circumstances the effects are mostly 
dominated by momentum considerations or not and to what extent initial values 
and the variability of the incident flow exert control on the dynamics. Many of 
the applications mentioned have complex interactions which are of multi-body 
type, with multiple impacts and rebounds of transported bodies, as well as of 
fluid-body-interaction type. It may be that considering the dynamics of a single 
body on its own is not immediately relevant to some of the real- world situations 
but on the other hand the single-body configuration does create a clear and 
basic growth point. A hope or hypothesis in principle is that this would develop 
fully to the multi-body arrangements. 

Clearly the criteria for lift-off of the body from the surface are a central 
issue to be addressed as part of the general picking-up of the body and trans- 
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porting it. This is so for many applications indeed but includes especially the 
matters raised in recent work about dust motions, saltation and reptation or 
splashing of particles on or near the surface of Mars. In particular, the first ex- 
traterrestrial measurements of sand transport rates were given in |Bridges et al.| 
2012a|b , and were rather surprisingly estimated to be comparable to those on 



Earth, in spite of the Martian atmosphere being 100 times less dense than the 
terrestrial one. These high transport rates occur despite only rare occasions 
on which the surface wind stresses on Mars are observed to exceed a supposed 



critical threshold for sand transport: Zimbelman 2000 , Haberle et al. 2003 



It has been hypothesised that the relatively low drag and gravity effects could 
keep sand in motion for longer as it bounces and tumbles across the surface [Kok 
2010a|b 2012 , to which the present study is directly applicable. Equally, we 



can here address the interesting question raised by Sullivan et al. 2008| as to 
whether the aerodynamic effects studied herein contribute to the lifting of Mar- 
tian sands Wang 2012 . More general work on dust suspended in the Martian 



atmospheric boundary layer is also of interest here Taylor et al. 2007 , Davy 



et al. 2009 . We shall return to discuss the Martian matters in some detail near 
the end of the present contribution. The phenomena involved in the broad area 
are dependent on many parameters for sure but are also dynamic, evolutionary 
and initial-value dependent to one degree or another. The present model which 
is depicted in figure [Ta| below and the ensuing working are based on consider- 
ing theoretically phenomena which are indeed dominated by momentum and 
pressure forces and hence centre on essentially inviscid unsteady fluid-body in- 



teractions of the type studied recently by Smith and Wilson 2011 , Hicks and 



Smith 2011 , Smith and Ellis 2010 in various contexts. These last papers 
find that due largely to the actions of added mass in the interactions induced 
between body and fluid interesting real-world phenomena such as touchdown, 
lift-off, clashing, skimming and rebounds emerge as part of the theory. Rela- 
tively few ingredients are necessary but the detailed evolution does tend to play 
a substantial role. 

The fluid is taken here to be Newtonian and incompressible with uniform 
density pjj say, where the subscript D refers to a dimensional quantity. The 
generally unsteady motion of both the fluid and the immersed body is assumed 
to be two-dimensional as a starting point for the theory even though it must 
be accepted that this restriction leaves out the distinct possibility of fluid skirt- 
ing around a contact point in the third spatial dimension. The representative 
Reynolds number Re based on incident flow speed and a typical body length 
varies from application to application above but the typical Re value of interest 
is quite large and so as a first approximation an inviscid separation-free theory 
is applied, in keeping with the overwhelming nature of the momentum described 
previously for these fluid-body interactions. Neglect of viscous effects seems ac- 
ceptable since most of the fluid-body interactions in reality are of the turbulent 
kind rather than the more sensitive laminar type. Thin bodies are also of inter- 
est in their own right because they can lead to fulsome analytical descriptions 
for many configurations such as in the previous paragraph but further they shed 
some light on interactions for thicker bodies such as smooth bluff cylinders or 
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spheres wherever thin layers of fluid lie between the body and the supporting 
solid surface. 

The paper itself is structured as follows. Section [2] considers the model in 
detail including particularly the generation of fluid-body interactions, followed 
by section [3] which turns to an analysis for relatively small times and investi- 
gates the possibility of an instant lift-off of the body from the solid surface when 
the oncoming fluid stream is instigated. Numerical studies for order-one times 
are presented in section [4] for cases where lift-off either is absent or possibly is 
delayed, with fluid-flow effects being suppressed. Sections [5j [6] then include nu- 
merically and analytically the added influences of fluid flow and the mechanism 
and criteria for lift-off to take place, respectively. Certain parameter effects are 
found to exert an important influence at that stage. There are indeed several 
parameters to take into account such as the scaled mass, moment of inertia, 
gravity, along with the scaled initial conditions. The final section [7] provides a 
further discussion. 



2 The fluid-body interaction 

The interaction between the body and the surrounding fluid is described as 
follows, with non-dimensional variables being used throughout. Cartesian ve- 
locity components {u, v) and pressure p are based on the horizontal incident 
fluid speed Ud and poU^, respectively, while the corresponding Cartesian co- 
ordinates (x,y) and time t are based on the characteristic body length Lp and 
travel time Ld/Ud in turn. The oncoming flow is thus a uniform stream with 
u = 1, v = 0. The body is thin, being of typical scale 0(h), h -C 1, in y and 
length unity in x, and is taken to be of smooth shape. In the configuration 
diagram of figure [Ta[ the stream is depicted as moving fluid from left to right. 
The surface, or wall, upon which the body rests is horizontal. The body is in 
contact with the wall, at least in the first instance, at a single contact point as 
shown. An inviscid or quasi-inviscid theory is applied here, with separation-free 
unsteady flow assumed in the two thin fluid-filled gaps labelled I, II on either 
side of the contact point, whose position varies with time. 

In the left-hand gap I ahead of the unknown contact point at x = Xo(t), we 
have in non-dimensional terms 

H t + (uH) x = , (la) 
u t + uu x = -p x , (lb) 

P+\ u2 = \ atx = 0, (lc) 
u = x' (t) at x = xo(t) . (Id) 

Subscripts t, x denote partial derivatives, whereas the prime denotes an ordi- 
nary derivative with respect to the relevant variable, in this case time t. The 



condition ( lc I allows for a jump across the leading-edge Euler zone, with the 



leading edge itself being at x = and with inflow being supposed locally. The 
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Figure 1: (a) A diagrammatic sketch of the flow structure showing the main two 
thin layers I, II under the rocking body, the fixed centre of mass (CoM) at x = x c , 
the contact point x = Xo(t) at the current time t, and the oncoming stream 
of fluid. Here in Cartesian coordinates, x is horizontal and y is vertically up- 
wards, (b) The four representative body under-surface shapes considered herein. 
The sinusoidal body has F(x) = sin(7ra;); the elliptical has F(x) = \J x — x 2 ; 
the smooth-cornered has F(x) — ~^X a {x) — ^Xb(x) + yJX a (Q) + y/Xb(x), 
where X a (x) — (x — ai) 2 + af, Xb(x) — (x — bi) 2 + b?,, and (01,02,61,62) = 
(0.2,0.1,0.8,0.1); and the constant-curvature has F(x) — x(l ~ x). 




attachment property (Id), which is discussed in some detail later, is associated 
with smoothness of the local flow solution close to the contact point. 

In the right-hand gap II downstream of the contact point we have similarly 
the governing equations 

H t + (uH) x = , (2a) 

u t + uu x = -p x , (2b) 

p = at x = 1 , (2c) 

u = x (t) at x = Xo(t) , (2d) 



differing from ([!]) only in (2c I. This condition at the effective trailing edge 
x = 1 in II is appropriate to the external-flow configuration, with atmospheric 
pressure being taken as zero for convenience. The local flow in the thin gap is 
supposed to be outward. The condition also applies, however, to internal flows 
when the gap on one side of the body (in this case, the lower side) is small 



relative to that on the other side, as in Smith and Ellis 2010 . The reason for 



the requirement ( 2c I is that under the present assumption of a long thin body 
the pressure varies typically by only a small amount of order h throughout the 
external flow compared with its characteristic 0(1) variation within the two 
gaps as shown in |1J,-d, [2^,-d). The same requirement results if the part of the 
body exposed to external flow is bluff and preserves the separation-free motion. 

Coupled with the fluid-flow equations above are the body-motion equations, 
namely 

Mh"(t) = / p(x, t) dx + N(t) - Mg+ , (3) 
Jo 

18" (t) = [ (x- x c )p(x, t) dx + (so - x c )N(t) . (4) 
Jo 

Here x = x c is the prescribed x-location of the centre of mass and h(t) its 
unknown vertical y- location, while 9{t) is the unknown angle the body chord 
line makes with the horizontal. The factors M, / are scaled effects of mass and 
moment of inertia, while g + is the scaled acceleration due to gravity, and N(t) 
is the unknown scaled normal reaction force acting on the body due to contact 
with the wall. The dimensional mass and moment of inertia are p^L^M/h, 
PdL a d 1 /h respectively, the dimensional gravity is hU^g + / Lo (the Froude num- 
ber is thus l/(hg + ) while the Richardson number is hg + ) and the dimensional 
normal reaction force is pdU^LjjN . In formulating the body-motion equations 
( 3p | we have supposed that N(t) is positive; if it should ever turn out that N 



becomes negative then the body will be taken to be no longer in contact with 
the wall. The marginal case N = remains moot. The integral contributions in 
( 3p | are dominated by the gap pressures for the assumed thin body but if the 



external part of the body is bluff instead then the latter part also contributes. 
An assumption of integrability also needs to be mentioned with regard to in- 
tegration through the xq contact point in (3]4). In addition, the unknown gap 
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shapes upstream and downstream of the contact point are given by 

H(x,i) = h(t) + (x - x c )9(t) -F(x) , (5) 

where F(x) is the prescribed smooth shape of the under-surface of the body. 
The four main shapes considered in this paper are shown in figure |lb| Both h 
and 9 are functions of t which are unknown in advance. Clearly at the contact 
point the constraints 

3H 

H = ,— = a,tx = x (t) (6) 
ox 

hold for the smooth shapes considered herein. 

Our task in general is to solve the nonlinear system ([l])~([6]) for u,p, h,9 in 
effect. Of interest first is the behaviour just after the motion starts at time 
t = 0, say. An analysis of the fluid-body interaction then, as presented in the 
following section, is found to complement the subsequent numerical work as well 
as to yield helpful results concerning the physical understanding of the nonlinear 
interaction present. 



3 Small-time properties 

The body is supposed to be positioned initially with its contact point with the 
wall being at some station xq — A say when the stream is suddenly switched 
on at time zero. This is equivalent to abrupt application of a fluid flow or a 
significant change in the fluid flow. The body motion itself is assumed to start 
from rest. A first guess for the response at small positive times is that the 
constant-pressure form where p = 1/2 throughout layer I, with p = in II, 
might work as an initial condition with zero initial flow in each layer, since then 
(lbc), (2b c) are all satisfied. The guess is then modified, however, by ( 3p l 



requiring the variation in h, 9 to be of order t 2 , implying a variation of order t 2 
in the gap thickness H from ([5]), which leads on to u being of order t via ( la 



and hence the pressure variation is of order unity. As a result, we expect the 
small-time expansions to take the form 

p = p (x) + tpi(x) + t 2 p 2 (x) + . . . , (7a) 

u = tui (x) + t 2 u 2 (x) + . . . (7b) 

in layer I, with the coefficients being functions of x to be determined, and 
similarly 

P = P Q (x)+tP 1 (x)+t 2 P 2 (x) + ... , (8a) 

U = tU 1 (x) + t 2 U 2 (x) + ... (8b) 
within layer II. The corresponding body movement has 

h(t) = h a +t 2 h 2 + ... , (9a) 

9{t) = 6 + 1 2 9 2 + . . . , (9b) 

x (t) = A + tB + t 2 C + ... , (9c) 
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where < A < 1, B, C are constants, but inspection of the requirement ( Id I for 
attachment soon establishes that B = 0. Here the initial state satisfies 



h - F{A) + {A- x c )6 = , (10a) 

F'(A) = o , (10b) 

by virtue of the contact condition ([6]). 

Substitution into (|6j now shows at order t 2 that the shape effects h 2 ,9 2l C 
must be connected by the two relations 

h 2 + (A - x c )6 2 = , (11a) 

CF"(A) = 9 2 , (lib) 



if the smooth body remains on the wall. An interpretation of ( 11a I can be made 
in terms of a combined upward and rotational movement of the body at small 
times. 

Meanwhile, the body motion in ( 3|4 | requires at leading order the balances 

2Mh 2 = ,h+N - Mg+ , (12a) 
2I6 2 = J 2 + {A- x c )N (12b) 

to hold, where No is the leading O(l) contribution to the reaction force N, and 
Ji , J 2 are defined as 

J\ = I Po(x) dx + / Pq(x) dx , (13a) 



J 2 = I (x — x c )po(x) dx + / (x — x c )Po(x) dx . (13b) 

JO J A 

Thus by elimination of iVo a relation between h 2 , 9 2 is inferred, namely 

a 2 h 2 + (i 2 6 2 = 72 , (14a) 

with 

a 2 = -2M(A - x c ) , (14b) 

ft = 21 , (14c) 

72 = J 2 + (A-x c )(Mg+ - ,h) . (14d) 

The coefficients a 2 , /3 2 are known constants, whereas j 2 clearly depends through 
J%, J 2 on integral properties of the unknown pressure coefficients po(x), Pq(x): 
this represents an influence of added mass. We must move on to determine 
those pressure coefficients in terms of h 2 , or 6 2 , or both. 

In the fluid flow in layer I, from (la) at the order of t we obtain an equation 
controlling u\{x), given H (x) = h + (x — x c )9 — F(x) as the known initial 
gap shape. Integration in x then leads to 

Ul (x)Ho(x) - -2 (h 2 + ( £±i - x c ) 9 2 ) (x-A). (15) 
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A requirement of finiteness for u\ at the original contact position x = A is 
imposed in order to keep the associated pressure coefficient finite at contact in 
line with the fluid-body interaction structure. In fact, (lb) at leading order 
indicates that the induced pressure gradient remains finite at x = A— since 
h 2 + (x — x c )8 2 is then of 0(x — A) from ( 11a I, while the gap shape H (x) is 
then of 0(x — A) 2 in view of (10 1. The results now yield 



ui(x) = - 
p (x) = 9 2 



2{x~Af 
H (x) ' 
x (x - A) 2 



dx+p (A-) 



(16a) 
(16b) 



H {x) 

where it is noted in particular that the integral is convergent for all x. Further- 
more, the leading-edge condition (lc) becomes po(0) = 1/2 here, which relates 
Po{A) to 2 and leaves us with 



where K\{x) 



x (x - A) 2 



dx 



(17a) 
(17b) 



H (x) 

for < x < A. Exactly the same approach applies in layer II except for Po 
replacing p an d the trailing-edge condition (2c I replacing (lc). Therefore we 
find 



where K 2 (x) 



p (x) = e 2 K 2 , 

x (x - A) 2 



H (x) 



dx 



for A < x < 1. The functions K\,K 2 are specified functions of x. 
Hence returning to ( 13|14 ) we have now 

a 2 h 2 + (3 2 8 2 = s , 

with j3 2 — (3 2 — r and the constants r, s are 



(18a) 
(18b) 



(19) 



r 



(x — x c )Kx(x) dx + / (x — x c )K 2 (x) dx 



-(A- x c ) 
A(A - 2x c ) 



f / K x {x) dx+ t K 2 {x) dx 

\J0 J A j 



+ (^4 



(20a) 
(20b) 



The main system to be tackled is therefore (11a), (19) for h 2 ,9 2 . Afterwards 
C is determined by (lib), and we can in addition use (12) to check on whether 
the primary reaction-force contribution Nq is positive or negative. It is note- 
worthy here that the contact-attachment requirement ( Id ) at the present level 
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of approximation reprodcues the smooth shape condition (lib I, since the cur- 
vature effects Hq and — F" are identical. Moreover, the solutions for po,Po in 
( 17|18 ) indicate that there is an O(l) pressure jump produced across the contact 
position. 

We show the prime features of the small-time solutions in figure [2] for the 
representative body shapes shown in figure |lb| Some conditions are found to 
yield the reaction-force contribution No being negative, which corresponds to 
the occurrence of an immediate lift-off. For an (under-) body shape which is 
sinusoidal, as in F(x) = sin(7rx) between x zero and x unity, figure 2a shows 
the results for the scaled leading-order vertical reaction force Nq versus the 
scaled body mass M as the parameter A governing the initial location is varied. 
Here the scaled moment of inertia I is taken to be equal to the scaled mass 
M for dcfinitcness; keeping the ratio I/M fixed is meaningful in the sense that 
the body shape remains fixed but the dimensional density and relative gap 
width can still be varied for example. For the horizontally symmetrical shapes 
investigated in this paper a value A greater than 1/2 means the body's leading 
edge is raised above the trailing-edge height. The results indicate that for all 
the values A studied there is a finite range of values M(= 7) for which Nq 
becomes negative, implying that for such values the body immediately lifts off 
from the wall. On the other hand, as M becomes relatively large, the value 
of Nq always becomes positive and increases linearly with M, in keeping with 
asymptotic behaviour concerning the loss of fluid-flow effects then; the interval 
within which fluid-flow effects do matter is actually quite small, in this instance 
being confined to M less than 0.1 roughly. There is in any case a critical value 
of scaled mass M, corresponding dimensionally to the existence of a critical 
fluid speed for lift-off for a given dimensional mass. Figure |2b| is then for an 
elliptical shape such that F(x) = \fx — x 2 for x in (0,1), leaving the body 
slope and curvature singular at the edge points, whereas the sinusoidal body 
has finite slope and zero curvature at the leading and trailing edges x = 0,1. 
Here almost all values of A studied produce a range of values of M in which 
lift-off can be inferred because Nq < 0. The case A = 0.8 is the only exception 
among those presented. Similarly, figure [2c] which is for a smooth-cornered 
shape, shows ranges of negative Nq and hence lift-off in almost all cases. Here 
F(x) = —y / X a (x) — y/Xb(x) + \/X a (0) + \/Xb(x), where the spatially- varying 
contributions are given by X a (x) = (x — a\) 2 + a 2 , X&(x) = (x — bi) 2 + b\, and 
the constants are (ai,a2,b 1 ,b 2 ) = (0.2,0.1,0.8,0.1). With comparatively large 
M(= I) in figures 2b (c) the asymptotes are linear as in figure 2a By contrast, 
with comparatively small or in effect zero M associated with dominant fluid 



influences, an explicit form of ( 11a ) ,( 19 1 can be found: this predicts that for 



a shape of constant curvature, for example, the critical of A is 2/3 at zero M, 
a critical value which is broadly in line with the results in figures 2a -(c). The 
critical value is also confirmed by the results of figure [2d] which are for the shape 
F(x) = x(l — x) of constant curvature. 

Figures 2d 2e indicate respectively the rather sensitive dependence of lift- 
off or its absence on the body shape and on the gravity factor. Concerning the 
latter, although most results given in this study take the value of g + as 10, we 
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(a) 



(b) 




(c) 




(e) 




Figure 2: Numerical solutions for the small-time system (llab),(19). Each 
figure shows for different body shapes F(x) the scaled leading-order vertical 
reaction force Nq versus the scaled body mass M for various values of the 
parameter A governing the initial location. Here, / = M for dcfiniteness. (a) 
Sinusoidal body F(x) — sin(7rx), x € (0,1); g + = 10. (b) Elliptical body 
F(x) = \Jx — x 2 , x G (0,1); g + = 10. (c) Smooth-cornered body F(x) — 
-y/xjx) - ^X b (x) + ^X a (0) + y/X b (x), where X a (x) = (x - ai f + 



2' 



X b (x) = (x - h) 2 + b 2 , and (ai, a 2 , h,b 2 ) = (0.2, 0.1, 0.8, 0.1), x G (0, 1); g+ = 
10. (d) Constant curvature body F(x) — x(l — x), x G (0,1); g + = 10. (e) 
The influence of the (inverse) Froude number on lift-off. Here, the body shape 
is the constant curvature case of figure 2d with A = 0.7, but g + is varied from 



10 down to zero. (/) As for figure 2e but over a wider range of values of M, I. 
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may begin to explore now the influences of the (inverse) Froude number. Figure 



2e has the value of A remaining at 0.7 but g + is varied from 10 down to zero. The 
range of conditions that lead to immediate lift-off is increased considerably with 
such decreasing of g + , a matter which is pursued further in detail in Appendix 
[X]and leads into the investigations in later sections on finite-time responses. 

4 The behaviour at 0(1) times with negligible 
fluid effects 

Without significant effects from the fluid flow, the controlling balances in |l])-(|6| 
over the time scale of order unity reduce to 

Mh"{t) = N(t) - Mg + , (21) 

ie"(t) = (x - x c )N(t) , (22) 

H(x, t) = h{t) + (x- x c )6{t) - F(x) , (23) 
OH 

H = , ^=0 at x = x (t) . (24) 
ox 

Here the unknowns are N,xq,H,9 as functions of time t only. The coupled 

system above corresponds to the relative mass and moment of inertia being 



large. An assumption of rolling is also noted. Elimination of N in (21), (22 1 
leads to the equation 

16" (t) = M{x - x c ){h" + g+) , (25) 



which couples with (23), (24} for (x ,h,9)(t) 



Applying ( 24 ) on the other hand yields the following successive relationships 



between h, 9, and Xq involving the body shape function F(x), 

h{t) = F(x ) - {x - x c )F'(x ) , (26a) 

h'(t) = -(x -x c )F"{xo)x' , (26b) 

h"(t) = -F"(x )x$ - (x - x c )(F'"(x )x'* + F"{x Q )x>>) , (26c) 



9 = F'(xo) , (27a) 
?' = F"(x Q )x / , (27b) 

'" = F"'(x )x' 2 +F"(xoK ■ (27c) 



Hence from substitution into ([25j) we obtain a nonlinear governing equation for 
xo(t) alone, which is 

ax'i + /3x£ = zg + , (28a) 
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(a) (b) 

Figure 3: Numerical solutions for the case of order unity time and negligible 
fluid effects ( 28a b) subject to varying initial location of the contact point, while 
initial contact-point velocity is fixed at zero. In all cases, g + — 10, x c = 0.5. 
Figures (o) and (6) show the time evolution of the contact position xq for varying 
initial conditions and body shape: (a) is for a sinusoidal body F{x) — sin(7rx), 
x £ (0, 1), while (b) is for an elliptical body F(x) — \J x — x 2 , x € (0, 1). 



where 

a =(jf + z2 )f"> P= -/'" + *(/" + */"') , z = (x -x c ), f = F(x ) 

(28b) 



The fluid-free cases in (28i,b) can be integrated once to the energy form, in- 



volving a constant of integration which is fixed by the initial conditions, 



*J 9 <8c) 

for any effective body shape f(z). Further integration to obtain z(t) and hence 
xo(t) explicitly is shape-dependent, although there are interesting limiting sit- 
uations such as for the nearly-cornered shapes addressed elsewhere. It can be 
shown also that the reaction force N remains one-signed in all these fluid-free 
cases. 

Numerical solutions of ( |28a[ b), obtained using Matlab's ® ode45 solver, 
checked by an independent solver, and subject to varying initial conditions, are 
presented in figures [3j [4] for two of the earlier prescribed shapes of interest: 
the sinusoidal shape F(x) = sin(7rx) and the elliptical shape F(x) = sfx — x 1 . 
In both cases, the value of g + was taken to be 10, while the ^-position of 
the centre of mass was fixed as x c = 0.5. The sinusoidal shape is shown in 



figure 3a for initial conditions in which the contact-point velocity £o = x 'o(fy 
is kept at zero, and the initial contact location xq(0) is varied over the values 
shown in the legend. In all cases except that in which xo(Q) = x c = 0.5, the 
solution evolves quite soon into an apparently periodic state in which the body 
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rocks. The cases xq(Q) — 0.4 and xq(0) = 0.6 remain near equilibrium as t 
increases with sinusoidal-like oscillations being observed about xq — 0.5 in a 
gentle rocking motion of the body. Values of Xo(0) further away from x c , on the 
other hand, provoke rapid rocking as these cases come closer to failure or lift-off 
in the sense that they yield comparatively rapid behavioural changes whenever 
the contact point Xo approaches either of the end points. However, the second 
shape, namely the ellipse, is found to be associated with comparatively slow 
responses whenever the contact point xq approaches either of the end points, 
as indicated in figure |3b| Gently-rocking sinusoidal-like oscillations about the 
equilibrium point at 0.5 again seem implied for initial conditions Xq(0) = 0.4 
and Xo(0) = 0.6, whereas values of xq(0) further away from x c lead to extreme 
rocking with relatively long periods in which xq remains close to zero or unity, 
accompanied by rather rapid increases or decreases of Xq in between. 

Analytical solutions are also possible for certain special shapes but it seems 
more useful here to examine the limiting responses indicated by the results in 



figures 3a (b) when the contact point is near an end point. In the locally near- 



flat scenario of figure 3a the limiting response occurs comparatively rapidly 
around some time t = to, such that time t = to + eT, say, where the parameter 
e is small and T is strictly of order unity, and x = eX to leading order near the 



leading edge. In normalised form, the governing equation (28a) or (28c) then 



leads to the local behaviour of the contact position being given by 

X 2 = \{T-Tof +T 2 , (29) 

7T 

with A > 0, T > 0, T being the constants A = g + x c {{I/M) + x 2 )- 1 , T = e^xo, 
while To is the a priori unknown scaled time at which xo reaches a minimum. 
In the elliptical scenario of figure |3b[ by contrast, the time scale is relatively 
long with t = e~iT and x = eX again to leading order for small e, so that from 



( 28a c) the scaled position X is found to be given by 



X = = . (30) 

(r-i-A(T-T )2) 2 

Here, A and Y are as before, To is in general a different value from the sinusoidal 
case, and T lies between T — y / 3/AT and T + ; close to the singular 

points there the growth of X matches to the behaviour in the relatively fast 
transition regions in which x becomes of order unity. Similar accounts apply 
near the trailing edge, where 1 — x — eX. The limiting responses just described 
are presented in figure |4j and they appear to reflect well the solution prop- 
erties found in figure [3] corresponding to rapid rocking and long-time rocking, 
respectively. 



5 Behaviour at 0(1) times with fluid effects 

Here we address the issue of what happens to the combined nonlinear fluid-body 
interaction if it persists over the time scale of order unity with significant effects 
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I — Asymptotic prediction with e =0 .3 1 | — Asymptotic prediction with e =0.3 1 
fc>u 




6 6.2 6.4 6.6 6.J 
t 



(b) 



Figure 4: Comparing asymptotic predictions with numerical results in the case 
of negligible flow effects, and with conditions as in figure [3] (a) The analytical 
limiting response ( 29 1 for the sinusoidal body in scaled (left) and unsealed (right) 
coordinates. The dots in the right-hand figure represent the numerical solutions. 
The time Tq at which xq attains its minimum was obtained from the numerical 



results in figure 3a (6) As for figure Hal but for the analytical limiting response 



(30) for the elliptical body in scaled (left) and unsealed (right) coordinates. 



Here, a single values of e is plotted to make the left-hand figure clearer, and 
since in the right-hand figure the asymptotic predictions from the three values 
of e used in figure [4a] would coincide. Once again, dots represent numerical 
solutions, and the time Tq at which xq attains its minimum was obtained from 



the numerical results, as shown in figure 3b 
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from the fluid flow being present. Instead of pi]), (22) we have now 



Mh"(t) = i x + N(t) - Mg+ , 
IB" (i) = i 2 + (x - x c )N(t) , 



(31) 
(32) 



with being the integrals from fluid-flow pressure effects in @,@ respec- 
tively, given explicitly in ( 34b I below. Elimination of N here leads to the equa- 
tion 

16" [t) =12- (x - x c )h + M(x - x c )(h" + .g+) , (33) 



rather than (25). Further, the successive relationships (26c), (27) between h,0 



and xo involving the function F still hold, and so the nonlinear governing equa- 
tion for Xo(t) now is 



ax'o + (3x'q = zg + + 



i 2 - Zll 



(34a) 



in this full fluid-body interaction, where, to clarify, 



1 1 



p[x,t) dx 



'2 



(x — x c )p(x, t) dx 



(34b) 



are the flow-pressure contributions. 



Numerical solutions of |34| were derived by adding the fluid effects of ( 34b ) 



(the added mass, as it were) in lagged style into (34a), then solving that for 



an updated xo, feeding this latest Xq into the fluid-flow calculation applied 



upstream and downstream of the xo-station to determine u from (la), (2a) and 



hence p from (lb c),(2b c), allowing the latest p to be fed into (34b), and so on, 
iterating per time step. An interpolation was employed in the flow calculation 
to handle the movement of the contact point x = xq with time, making use of 



the requirements (Id), (2d) and enabling the quasi-mass flux to be evaluated. 



The typical values of the uniform time step St and spatial step 5x used in the 
computations were 0.001 and 0.005 in turn, and the effects on the solutions of 
varying the steps were tested thoroughly. The major results are presented in 
figures [sHT] for various prescribed shapes F and varied other conditions. 

Figure[5] is for the sinusoidal shape, and M = I = 0.125. Figures [5^,,b, 
in which the initial conditions are of Xq being 0.33 while the contact velocity 
£o is zero, show details of the evolving pressures and velocities in the fluid- 
filled gaps over a scaled time range of 0-10. Thus figure [5a] has the results 
for pressures p\ , pi ahead of and behind the contact point plotted against x at 



integer times, and figure 5b gives the respective fluid velocities u\,U2 versus 
x . It is interesting that the pressure curves upstream and downstream of 
contact are almost straight with only a slight curvature being apparent. Also the 
vertical lines in figure [5a] mark the clear jumps in pressure at the moving contact 
position, in keeping with the earlier remarks. The fluid velocities ui,U2 in 
contrast seen in figure 5b are continuous through the contact position x = Xo(t) 



as anticipated previously. Further, the velocities are all positive at early times 
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(c) 



Figure 5: Numerical solutions for the behaviour at 0(1) times with fluid effects. 
Here, the body shape F is sinusoidal, and M = I = 0.125. (a) Pressures pi,P2 
respectively ahead of and behind the contact point for the initial conditions 
(^O)Co) = (0.33,0). The discontinuity across the contact point is clearly visible. 
(b) Velocities u\,U2 respectively ahead of and behind the contact point for the 
initial conditions (xo>£o) = (0.33,0). The continuity and relative smoothness 
at the contact point are noted, (c) Contact location xq, its velocity £o, its 
local shape response F(xq), and the reaction force N versus t, for the initial 
conditions (xo,£o) = (0.33,0) (solid lines) and (0.28,0) (dashed). 





(a) 



(b) 
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(d) 

Figure 6: Numerical solutions for the behaviour at 0(1) times with fluid effects. 
Here, the body shape F is elliptical, (a) Here, for the initial condition (xq, £o) = 
(0.25,0) and for M = I = 0.08, are shown the contact location, its velocity, its 
local shape response, and the reaction force plotted against scaled time, (b) 
As for (a) but with M = I = 0.05. (c) Results for Xq,N for three different 
initial contact velocities £o equal to —0.2, 0, 0.2, with the same settings as in 
(b). Responses for the first value calculated with three different time steps are 
also shown (labelled "coarse grid", "medium grid", and "fine grid"), (d) The 
evolution of F(xq) and 9 with time are shown for the same conditions as in 
figure |6a 
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Figure 7: Numerical solutions for the behaviour at 0(1) times with fluid effects. 
Here, the body shape F is the smooth-cornered body of figure [2cJ with M — 
I = 0.1 and initial conditions (xo,£o) = (0.25,0). Beyond the lift-off time when 
N first becomes negative, numerical results are shown by dashed lines since 
strictly they are no longer physically meaningful even though interesting. 



but eventually negative values are encountered (smoothly) at later times as the 
front of the body rocks downwards squashing the fluid there in a sense; the 
inferred flow closer to the leading edge is then clockwise around the leading 
edge. Figure [5c] presents the numerical results for the evolution of the contact 
location xq, its velocity £o, its local shape response F(xo), and the reaction force 
N versus t, both for the initial condition (£Co,£o) of (0.33,0) as above, and for 
(0.28,0) for comparison. The latter case shows a phenomenon of rapid rocking 
appearing near the end points which is not dissimilar to that encountered in 
no-fluid scenarios subject to different initial conditions. 

In figure [6] the results for the elliptical shape are presented. Figure [6a] with 
an initial condition (xo,£o) °f (0.25,0) is for M = I = 0.08 and shows the 
contact location, its velocity, its local shape response, and the reaction force 
plotted against scaled time. The findings are fairly similar to those for the 
no- fluid case in figure [3] but clearly the fluid effect is no longer negligible here, 
although an apparently periodic rocking behaviour emerges nonetheless. Figure 
[Gb| has the same initial condition but now M = I = 0.05, implying increased 
fluid-flow effect. Here the contact position xq is found to increase monotonically 
with time, and tends to approach the trailing edge, but during that process the 
reaction force N becomes negative and so lift-off of the body is indicated. This 
arises at a finite time t which in the case of figure 6b is at t = 1.921 and the final 
contact at that time is at an xo value of 0.9096. Figure [6c] reinforces the point 
and further examines the effects of the initial conditions by showing the results 
for xo, N for three different initial contact velocities £o equal to —0.2, 0, 0.2, with 
the same settings as in 6b and with the responses for the first value also being 
calculated with three different time steps to highlight the accuracy involved. In 
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Figure 8: The influence of the gravity parameter g + (value in parentheses in 
legend) on the behaviour of the constant-curvature body of figures 2d 2lj Here, 
M = I = 0.1, and the initial conditions are (xo,£o) = (0.7,0). 



every case the reaction force becomes negative within a finite scaled time. The 
evolution of the angle 9 under the conditions of figure [6a] is presented in figure 
|6d| and also supports the conclusions on nonlinear periodic rocking behaviour if 
lift-off is absent, as well as agreeing with the small-time analysis of section [3j 
In figure 6d the successive maxima in F versus time are seen to correspond to 
9 changing sign, whereas the minima in F correspond to successive extrema in 
9 because of the relationship between F and the contact location Xq. 

The smooth-cornered shape of figure [2c] is the subject of figure [7] with fluid- 
flow influences now being present. The trend observed is akin to that described 
in the previous paragraph, including the monotonic response of the contact 
location prior to the eventual lift-off inferred from the reaction force N becoming 
negative at a finite time. Numerical results generated from the lift-off time are 
shown by dashed lines to indicate that they are not strictly physically meaningful 
even though interesting. 

The influence of the parameter g + was examined first in figures 2e and [2f] and 
is re-examined in figure [8] Reducing g + from its usual value of 10 to the values 
5 and then 3 is found to result in a change from periodic rocking behaviour to 
finite-time lift-off as N becomes negative as in figure [8] doing so earlier for the 
value 3 (at t = 0.817) than for the value 5 (at t = 3.576). The contact position 
moves monotonically towards the trailing edge for the two reduced g + values, 
and the associated 9 angles decrease monotonically, until lift-off is encountered. 
The lift-off is perhaps the most intriguing phenomenon to explore next. 



20 



6 Lift-off of the body 



The evidence from the results in figures [6j|8] indicates that the interactive solu- 
tion is regular at the onset of lift-off of the body, with (to repeat) lift-off being 
identified with the reaction force becoming zero after being positive beforehand. 
Thus the expression for the flow and body-movement behaviour at the lift-off 
time <lo say has 



oo 

(h,e,u,p,xo,N) = y(h,e,u,p,x ,N)^(t L o-t) m . (35) 



m=0 



The coefficients etc on the right-hand side for m ^ are independent of 



t but their values depend on the evolution from the initial time, while (35 1 holds 
for scaled time t tending to £lo~ > the lift-off time being positive. It is noticed 
from figure [la] and so on that the flow velocity u and pressure p are in two parts 
I, II on either side of the contact point. For the pressure in particular this facet 
implies there is a moving discontinuity to be dealt with as highlighted by figure 



5a and hence a little delicacy near the terminal location x = xo(£lo) a ^ sman 
relative times t^o — t, but this does not affect the response ( [35] ) materially. 
Substitution of (351 into the fluid-body equations of section [2] or into the 



equivalent system ( 34a| b) can then be made and studied in detail. The principal 



point however is that in (34a b) the integrated flow-pressure contributions i\, %2 
act to counter the effect of scaled gravity on the right-hand side and thus take 
the interactive solution into a part of solution space which is not attainable in 
the no-fluid scenario of section [H The action which is due to added mass and 
to evolution in the fluid-body interaction allows the reaction force to become 
negative. The same global type of anti-gravity action is also responsible for the 
changes of sign in the reaction force at early times (section [3]) as shown in the 
results of figure [2j 



The results for the ellipse in ( 30 ) , figure 3b without fluid effects and subse- 
quently figures |6^a-c) with fluid effects provide a specific example of the above. 
Checking on the orders of magnitude of gap properties and flow behaviour lo- 
cally near the contact point whether close to the leading edge (x = eX) or 
to the trailing edge shows dh/dt and dO/dt to be of typical size e~3 which is 
large. So V is of the same typical size. The fluid velocity u is then of order 
dxo/dt, that is, eJ, making q also small of order e5 . The local pressure response 
p is therefore merely a small perturbation from the value 1/2 and hence has 



negligible influence on ix, ii in (34ab) as these remain overwhelmed by global 



contributions of order unity. Similar reasoning applies for the sinusoidal shape 



of ( 29 1 and figures pal p[a-c) for example 



The conclusion then is there is nothing dramatic about the onset of lift-off, 



in line with ( 35 1 . In contrast such a conclusion is not necessarily so once lift- 
off occurs for times t equal to £lo + : the latter is considered within the next 
section. As regards the work so far the more significant question is whether 
lift-off occurs or not and this is clearly reliant on the interaction parameters 
together with initial values. According to the model lift-off generally cannot 
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take place without fluid flow effects but it can with them. 

7 Comparisons and final comments 

The lift-off of a body from a fixed solid surface due to fluid motion or indeed 
just the washing or rocking of such a surface-mounted body clearly depends on a 
significant number of parameters. The total parameter space should also include 
the influences of evolution (the initial-value problem due to wash starting up or 
wind changing for example) and body shape because of their importance. The 
present investigation it is hoped helps to shed light on this parameter space by 
means of the specific studies in sections [3}j5] on small-time responses, zero-fluid 
evolution and with-fluid evolution in turn, given the model set up in section [2] 
and the account of lift-off in section [6] In addition concerning scaled gravity 
effects in particular at low Richardson number two modes of lift-off can now 
be identified. One occurs for enhanced values of the scaled mass and so has 
the gravitational force in balance with the mass-acceleration contribution on 
account of a reduction in the local variation in gap width. Typically this mode 
corresponds in dimensional terms to increased incident fluid speed (including 
threshold speed) or increased body mass. The second mode occurs for in effect 
zero scaled gravity in the sense that the gravitational force simply has negligible 
influence on the mass-acceleration balance compared with the flow forces and 
the normal reaction force. Such a second mode can be associated mainly with 
relatively reduced body or particle dimensions or negligible gravity in reality. 
The specific exact case addressed in appendix [A] backs up the presence of these 
two modes of lift-off. 

It is interesting to return now to consider the movement of sand or dust on 
the surface of Mars, a matter introduced in section [T] As a broad background 
observations have suggested that the flux of Martian dust or sand movement is 
comparable overall to that on Earth. There are, obviously or potentially, many 
physical factors at work here as is indicated also in the substantial fairly recent 
growth in studies: again see section [T] The present study suggests in fact there 
are wide areas of parameter space where lift-off can occur in practice. 

Two attributes which might be emphasised first are that the Martian grav- 
ity go is about 0.38 of the Earth value whereas the density of the Martian 
atmosphere po is only about 1/60 of Earth value. Predictions then from 
our results would be based on knowing that M — M D h/{p D L D ) and g + = 
gu^D I '(hUp) from the non-dimcnsionalisation, and estimating the body mass 
as Mjj = pDhodyL 2 D x (unit distance in Zd), where pDbody is the typical body 
density. In consequence 

M = P^Lh 2 . (36) 

PD 

Thus on Mars, taking the same size, shape and density of body (dust particle) 
as on Earth, one should tend to find M increasing by virtue of the po factor 
in M although readily mitigated by the gap width h reducing by a factor 1/4 
say, making M increase by about 4. Now consider the mode mentioned earlier 
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for lift-off assuming M is relatively large, that is, M ~ l/.9 + or M < l/g + - 
This becomes the requirement Up > (pDbody/pD)gDhL]j on the threshold wind 
speed Ud for lift-off to occur. In numerical terms the right-hand side here as 
far as Mars is concerned is approximately (60) x 0.4 x 1/4 x 1 relative to the 
Earth value in view of the estimates above, i.e. the Mars value is about 6 times 
the Earth value on the right-hand side. Hence the critical wind speed on Mars 
is predicted as about 2-3 times that on Earth. This last result is reasonably in 



line with the references, e.g. Wang 2012] ; further the current prediction of a 



square-root dependence of threshold speed on particle size is not out of keeping 
with figure 2 of that last paper. 

Second, though, is the need for caution, since the above prediction and in- 
deed earlier predictions should be qualified heavily due to the following aspects. 
As we have shown, lift-off can appear quite readily at low scaled mass or large 
scaled mass typically: see also figure [2eJ Moreover shape effects, initial con- 
ditions and gap width all matter considerably as shown in the current study. 
Many areas of the total parameter space allow lift-off to occur. Again in reality 
the 'typical' body shapes vary appreciably rather than being purely spherical or 
thin, while with many bodies or particles present there is a potentially compli- 
cated and subtle multi-body factor where for example a string of bodies lies on 
top of another string, creating a thin gap, and we should not forget the existence 
of a laminar or turbulent boundary layer in the incident fluid motion. Finally 
here, as shown in the previous section, lift-off itself is not a dramatic event 
in terms of the fluid-body interaction per se: it merely happens in a smooth 
manner, at least in the pre-lift-off stage (compare with the post-lift-off stage 
addressed below). This manner is quite in keeping with the stated importance 
of the initial conditions. Parameter studies of rcptation in terms of fluid-body 
interaction along lines similar to those here or in Smith and Wilson 2011 could 
be of much interest. 

Intriguing issues remain. The influence of incident shear for example in the 
depths of an oncoming atmospheric boundary layer has still to be considered. A 
classic treatment of this aspect is for a uniform shear. In addition the modelling 
so far has ignored viscous effects quite generally. Other shapes or configurations 
of body need to be examined such as cases with non-symmetry in the stream- 
wise direction or with irregular shapes or groups of bodies, in order to raise 
understanding for dust piles for instance. (Appendix [A] shows that the body 
curvature can act in effect to increase the scaled mass parameter and reduce the 
gravity force, which leads on again to another part of parameter space.) 

Properties arising after a lift-off are also intriguing. In cases where lift-off 
does occur (negating the energy-integral result of the no-flow case) the local 
scales soon after lift-off are expected to be similar to those in the current study, 
such as a relative time squared scaling in the vertical coordinate. This would 
suggest that a local analysis then has to deal with a substantial jump in pressure 
within which the local pressure is an order-one function of distance measured 
from the lift-off point but scaled with relative time. In some detail near a moving 
contact point prior to lift-off of a smooth body the local gap shape is typically 
of the form H ~ (j,(t)x 2 if the contact point is at x — Xo(t). Here fi(t) is an 
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O(l) function of time which is arbitrary in the sense of being determined by 
historical effects, and x denotes x — xq{£). Now the response in horizontal fluid 
velocity is of the form u ~ x' + rx where the rate r(t) is /// (3/x), implying that 
the pressure behaves as c 2 (t) + x'^x + (r' + r 2 )x 2 /2 in the vicinity of the contact 
point. However the c 2 pressure function is discontinuous across the contact 
position, as is clear from the results in figure [5a] for example and the working in 
earlier sections. This reinforces the suggestion of a pressure discontinuity which, 
immediately after lift-off, must be accommodated by the fluid flow through the 
newly opened gap. 

The study has been focussed throughout on unsteady nonlinear interactions 
in just two spatial dimensions. Moving the theory on to three-dimensional 
interactions is called for. The three-dimensional version of the post-lift-off con- 
figuration described in the previous paragraph for example would be of much 
interest. Many other extensions naturally suggest themselves: allowing body 
flexibility (with applications to red blood cell deformation and swimming mi- 
croorganisms amongst others); studying the effects of surface shape, curvature, 
and roughness; and including the effects of a thin layer of a second, viscous fluid 
between the surface and rocking body, to study the effects of surface tension 
there. 

Thanks are due to members of the Medical Modelling Group at UCL and the 
UK Icing Group for their interest and comments, and to Dr Miguel A. Moyers 
Gonzalez at UC for useful discussions. 



A The constant-curvature body at small times 



A constant curvature in the present setting corresponds to having the initial 
gap width 

H (x) = X(x - A) 2 , (37) 

where A is a positive constant proportional to the body curvature. Hence eval- 
uating the integrals in ( 1 7a)-( 20a| and in particular finding that r is (— A 2 + 
A - l/3)/(2A) leads to the result 



A 2 -A 



2A 



+ 2I + 2M(A- 



(38) 



for 9 2 in terms of s which is Mg + (A — 1 /2) — A 2 /A. The value of x c is taken as 
1/2 as a major example. Then the scaled reaction force No follows from (12b I 
in the form 

No = "-^f . (39) 



Here J 2 can be determined analytically from ( 13b) with ( 17a I, ( 18a I and so this 
yields the explicit solution 



No 
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A —A 
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(40) 
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Figure 9: The constant curvature body at small times, (a) Solutions of equation 
(40) for a range of values of g + 



9a 



but for a wider range of 
/. (c) The numerical roots of (40) obtained from figures [9a| |9b| 



(b) for figure 



values of M 

(circles) compared with the asymptotic prediction of equation ( 42 ) (solid curve) 



for No, with 62 given explicitly by (38). It can be shown by working with the 



quantities AM, 82/ X, g + /A that A may be normalised to unity without loss of 
generality. 

Plots of the scaled force N against the scaled mass M are presented in figures 
9a and 9b for varying values of the gravity contribution g + . In these plots the 
scaled moment of inertia / is again equated to M for reasons described earlier 
on and the curvature constant A is kept at unity while the starting location 
A is kept at 0.7 for comparison purposes. The results here agree with those 
presented in section [3j 

The lift-off case where N is zero is of much interest of course. This critical 
case is controlled by the relation 



2/jg+M 2 



A(A-l)(A-h) , (A*-A + kW 



{A- 



l)A 



^ + — T " + V " 2a'^ M+ '^sf ZZ = 

(41) 

between M, g + , A and A, from (40) with (38), with \i standing for the ratio 



I/M. When the gravity effect is comparatively small for instance the first two 



contributions are overwhelming when M is large and the relation (41) leads to 
the explicit form 



M . 



[2 M +(A-l)(A-i)]A 
4 W + 



(first root for g + small) 



(42) 



for the dependence of the critical M value on g + in particular. The numerator 
here can be positive or negative, depending on the ratio fj, and the contact 
position A, which indicates a sensitivity to the precise body shape. However for 
the values [i — 1 and A — 0.7 used in figures [9a| and [9b| the asymptotic trend 



( 42 ) is found to agree rather closely with the full numerical results in the figures 
as g + varies and indeed the asymptote even works well for g + equal to unity or 



greater, as indicated by figure 9c . The trend (42) for this special category of 



body shapes also hints at the trend for the general body shape, which again is 
found to have M proportional to the inverse of effective gravity. 

Similar considerations apply to the other root for relatively small gravity in 



(41), giving the explicit value 



M . 



(A- 



l)A 



8\[/iA 



A(A-l)(A-i), 



(second root for for g + small). 



(43) 



This is associated with having identically zero gravity in (41) and corresponds 



to one of the cases presented in figure [9j The potential reliance on A exceeding 
2/3 is again noted. 
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